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y  The  hvdrodvnamic  processes  resulUng  from  instantaneous^  beating  of  a  cylindrical  region  of 
gas  are  studied  by  means  of  numerical  integration  of  the  fluid  equations  Various  iniUal  over 
pressures  and  radial  profiles  of  initial  presaute  are  considered.  For  very  smooth  radial  profiles,  eg. 
the  Rennett  profile,  the  evolution  is  verv  simple:  the  gat  expands  radtallv.  a  single  outgoing  shock 
forms,  the  density  at  late  times  undershoots  the  value  at  pressure  balance  by  only  about  HPT  and 
then  increases  to  pressure  balance  tiotsly .  monotonicallv .  and  atfiabaticaily .  For  profiles  that  are 
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p  more  square  in  ahape.  or  are  aharpty  diaconlinuoua.  the  evolution  i*  much  mote  complicated.  A 
tiron|  rarefaction  typically  occur*,  of  bnef  duration  and  in  a  narrow  teflon  about  the  axis.  Durtnf 
the  rarefaction,  both  the  density  and  temperature  fall  welt  below  the  value*  at  pleasure  balance.  The 
rarefaction  is  abruptly  terminated  by  an  inward-propagalinf  shock.  followed  by  multiple  drnaity 
oartllattons.  which  can  form  (hocks  for  ttronf  overpressure*.  The  dependence  of  the  time  scale*  and 
strengths  of  the  rarefactions  and  shocks  on  the  initial  overpreaaure  and  profile  shape  are  calculated 
and  discussed 
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GAS  CHANNEL  FORMATION:  DEPENDENCE  OF  HYDRODYNAMIC 
PHENOMENA  ON  THE  INITIAL  RADIAL  PRESSURE  PROFILE 


I .  Introduction 

Several  laboratories  have  recently  studied  the  production  of  hot 

reduced  density  channels  in  air  and  other  gases  by  pulsed  directed 

1  2 

energy  sources  such  as  electron  beans ,  pulsed  lasers,  electrical 

3  1* 

discharges,'  and  exploding  vires.  Such  channels  persist  in  pressure 
balance  with  the  surrounding  gas  until  they  cool,  typically  for  milli¬ 
seconds  and  longer.  All  of  the  energy  sources  cited  above  typically 
have  pulse  durations  ( <<  1  ^isec)  short  compared  to  the  hydrodynamic 
response  tine  of  the  heated  gas.  Thus  the  hydrodynamic  response  de¬ 
couples  from  the  energy  deposition  process,  and  can  be  treated  as  a 
one- dimensional  initial  value  problem:  the  radial  motion  of  an  ini¬ 
tially  .r.ifcrm- density  ga3  with  a  given  cylir.drically  symmetric  initial 

temperature  profile.  Self-3imilar  blast  wave  solutions  to  this  problem, 

5 

in  closed  form,  are  veil-knovn,  but  only  for  very  large  initial  over¬ 
pressures.  In  the  range  of  greatest  interest,  overpressures  ranging 

o 

from  tvo  to  10  times  the  ambient  pressure,  numerical  solutions  are  re¬ 
quired  in  general. 

g 

Vie  have  previously  examined’  the  hydrodynamics  of  channel  formation 
for  the  special  case  of  a  Bennett  initial  radial  overpressure  profile, 

typical  of  a  uniform  gas  heated  by  a  self-pinched,  mcno-energetic 
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electron  bean.7  The  evolution  was  found  to  be  remarkably  simple  in 
that  case,  with  the  central  heated  region  expanding  adiabatlcally ,  and 
nearly  monotonlcally ,  until  it  reaches  pressure  balance  with  the  sur¬ 
rounding  cool  gas.  The  central  density  dips  only  about  10Z  (practi¬ 
cally  Independent  of  the  overpressure)  below  its  value  at  pressure 
balance  and  then  very  slowly  increases  until  it  returns  to  pressure 
balance.  No  density  oscillations  occur  during  this  process.  A  single 
shock  propagates  outward  radially,  heating  the  gas  it  overruns,  and 
thus  moderately  reducing  the  final  density  in  the  radial  wings  of  the 
channel  and  broadening  the  wings,  but  having  no  effect  in  the  channel 
center.  Formulas  for  the  final  density  profile  (including  shock  wave 
effects),  and  for  the  time  evolution,  were  derived  heurlst lcally ,  and 
were  found  to  agree  well  with  numerical  solutions  of  the  hydrodynamic 
problem. 

Heating  profiles  differing  significantly  from  Bennett  can  result 
from  several  of  the  pulsed  energy  sources  mentioned  above.  In  parti¬ 
cular,  energy  deposition  by  lasers  through  clean  air  breakdown  occurs 
only  when  the  laser  Intensity  exceeds  a  sharp  threshold;  thus  the 
central  region  of  the  pulse  can  heat  the  gas  strongly,  while  the  radial 
wings  do  not  deposit  any  energy,  leading  to  a  nearly  square  or  "top-hat" 

g 

radial  profile  of  gas  temperature.  In  this  report,  we  present  detailed 
numerical  studies  of  the  hydrodynamic  evolution  resulting  from  a  variety 
of  radial  heating  profiles,  ranging  from  Bennett  and  Gaussian  to  square. 

Q 

Our  results  confirm  the  qualitative  effects  studied  elsewhere  by 
approximate  numerical  methods,  but  are  more  extensive  and  numerically 
accurate.  We  find  that  heating  profiles  with  a  more  square  shape  than 
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the  Bennett  profile  typically  lead  to  multiple  oscillations  (ringing) 
of  the  channel  density.  If  the  overpressure  is  strong,  these  oscilla¬ 
tions  can  steepen  into  shocks,  and  one  or  more  ingoing  shocks  can  over¬ 
run  the  ch*-.nel  and  reflect  off  the  axis,  in  addition  to  the  initial 
outgoing  shock.  For  brief  periods  during  these  oscillations,  a  strong 
ir.vard-propagat ir.g  rarefaction  wave  drives  both  the  density  anu  the 
temperature  or.  the  channel  axis  to  veil  below  their  final  values  (when 
the  channel  reaches  pressure  balance  with  the  ambient  gas).  (For 
example,  for  a  square  profile  or.d  initial  overpressure  factor  23,  the 
minimum  channel  density  is  65  of  the  initial  density,  while  the  final 
density  i3  95  of  the  initial.  Thus,  transient  situations  occur  that 
are  particularly  favorable  for  applications  where  a  cool,  low  density 
channel  is  desired. 

In  lection  2,  ve  present  the  results  of  numerical  calculations 
over  a  broad  range  of  overpressures  and  radial  profiles,  emphasizing 
the  dependence  of  density  minima,  shock  characteristics,  and  time  scales 
on  these  parameters.  The  results  are  interpreter  in  fee.  3,  and  cur 
conclusions  are  summarized  in  Sec.  - . 

2 .  Numerical  Results 

Ve  have  used  a  or.e-dimer.sior.al  (radial)  fluid  code  to  calculate 

the  time  evolution  of  the  channel  hydrodynamics  for  a  variety  of 

heating  strengths  and  profiles.  The  code,  which  has  been  described 
9,10 

previously,  perfcrm.3  finite-difference  solutions  of  the  equations 
of  conservation  of  mass,  momentum,  and  energy,  using  an  Eulerian  flux- 
corrected  transport  algorithm,  which  accurately  resolves  shocks.  All 
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(4a) 


?  (7,  7  -  0)  -  1  +  Pq(1  +  7  n)"n 

and  "generalized  Causslan"  profiles, 

P  (7,  7  -  0)  -  1  ♦  P  exp  (-  7  ~n).  (4b) 

The  value  n  ■  2  corresponds  to  Bennett  and  Gaussian  profiles,  in  Eqs. 
(4a)  and  (4b)  respectively,  and  Increasing  values  of  n  correspond  to 
squarer  profiles;  for  n  -*  «*  ,  both  profiles  approach  step  functions 
(but  the  generalized  Gaussian  profiles  approach  this  limit  faster). 

The  generalized  Bennett  and  Gaussian  profiles  are  Illustrated  in  Figs. 
1  and  2  for  a  wide  range  of  n. 

For  purposes  of  comparison  with  our  previous  work^,  we  show  here 
in  Fig.  3  the  time  evolution  of  a  deposition  with  the  Bennett  profile 

with  P  ■  37, for  earlv  times.  The  density  at  the  center  decreases 

o 

monotonically  and  smoothly  (the  only  shock  forms  at  r  •  .8  and  runs 
outward).  Th.'s  result  is  typical  for  the  Bennett  distribution 
profile  where  no  ringing  of  the  central  channel  density  is  observed 
and  no  ingoing  shock  appears. 

The  evolution  becomes  increasingly  complex  as  one  considers 
initial  temperature  profiles  that  are  squarer  then  Bennett,  as  shown 
in  Fig.  4,  where  the  time  evolution  of  c  (?  •  0)  is  plotted  for  a 
series  of  generalized  Gaussian  profiles  with  increasing  values  of  the 
index  n,  all  with  PQ  •  37.  Even  for  a  Gaussian  profile  (n  ■  2),  there 
is  a  gentle  maximum  of  7  (7  •  0,  7)  at  7  •  1,  and  for  n  t  3,  thla 
maximum  clearly  corresponds  to  an  Inward-propagating  shock  striking 
the  axis.  The  strength  of  this  shock  Increases  for  increasing  values 
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of  n,  until  (surprisingly)  it  reaches  a  maximum  at  n  as  10,  with 
density  jump  by  a  factor  %  7.3,  and  weakens  slightly  to  density  jump 
%  4.6  for  still  squarer  profiles  with  n  i.  10.  (The  fact  that  the 
density  jump  exceeds  the  maximum  Rankine-Hugoniot  value,  6,  for 
y  •  7/5,  is  due  to  the  combined  effect  of  the  inward  shock  and  the 
subsequent  outward  shock  reflected  off  the  axis.) 

For  n  >  3,  the  initial  pressure  profile  is  sufficiently  square 
that  no  motion  occurs  at  ?  ■  0  until  a  well-defined  rarefaction  wave, 
traveling  at  the  local  sound  speed,  reaches  the  axis.  The  density 
then  falls  rapidly,  reaching  a  minimum  value  that  is  as  low  as  J  ■ 
0.06  for  Pq  l  10,  Just  before  the  Inward  shock  reaches  the  axis.  For 
purposes  of  comparison,  adiabatic  expansion  would  reduce  the  density 
only  to  o  •  0.2,  for  ■  10.  The  duration  of  the  density  minimum  is 
quite  short,  about  4t  %  0.1.  Also,  for  n  i  3,  a  weak  but  definite 
ringing  of  the  channel  occurs  subsequently,  the  amplitude  of  which 
increases  as  the  heating  profile  becomes  squarer.  Such  ringing  might 
possibly  have  some  Impact  on  channel  stability  and  cooling  rate. 

In  Fig.  5,  ?  (r  •  0,t)  is  plotted,  for  various  overpressures  P 

o 

for  an  essentially  square  radial  profile  (generalized  Gaussian  with 
n  ■  50).  We  see  that  an  initial  ingoing  shock  is  well  formed  over 
the  entire  range  of  pressure  Pq.  The  arrival  time  on  axis  is  a 
function  of  the  overpressure  Pq  but  appears  to  be  nearly  Independent 
of  the  squareness.  The  strength  of  the  shock  and  the  value  of  the 
density  minimum  for  the  preceding  rarefaction  is  strongly  dependent 
upon  the  squareness  but  only  weakly  dependent  on  P  .  In  addition  we 
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see  from  Fig.  5  as  Is  increased  a  second  density  maximum  arises 
which  steepens  Into  a  shock  for  Po  i  56. 

A  series  of  plots  In  Fig.  6  of  p(?  -  0,  t)  for  generalized  Bennett 
profiles  with  various  values  of  n  but  at  Pq  •  74  shows  phenomenology 
similar  to  that  seen  with  generalized  Caussian  profiles  in  Fig.  4.  We 
see  here  that  the  second  density  maximum  steepens  Into  a  shock  for 
n  i  4.5. 

In  Fig.  7,  we  show  a  series  of  snapshots  at  different  times  of  the 
Instantaneous  density  profile  p(r),  for  a  generalized  Caussian  profile 
with  n  -  50  and  P^  •  74.  For  the  latter  case,  we  note  particularly  the 
profile  ?(r)  during  the  transient  density  minimum.  Just  before  the  In¬ 
ward  shock  arrives  on  axis.  The  low  density  "hole"  on  axis  is  rather 
narrow,  extending  only  to  ?  ^  0.3,  and  its  duration  Is  only  about 
5?  *  0.1,  as  the  hole  is  abruptly  closed  by  the  incoming  shock.  The 
depth  of  this  minimum  depends  quite  strongly  upon  the  squareness  of  the 
profile  (the  index  n),  as  seen  in  Fig.  4,  but  is  quite  insensitive  to 
the  initial  overpressure  as  can  be  seen  from  Fig.  5.  Thus  a  deep 
density  channel  can  be  produced,  momentarily,  even  by  a  relatively  weak 
beam.  However,  this  economy  of  means  might  be  difficult  to  use  for 
practical  applications,  because  of  the  brevity  of  its  duration  and  the 
precise  timing  required. 

3.  Discussion 

We  have  shown  previously6  that  the  channel  density  profile  at  late 
times,  when  the  channel  reaches  pressure  balance  with  the  ambient  gas, 
is  defined  by  the  set  of  equations 
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fT  (7,  *t  *  »  3  [?  fr  Cr) ,  ~  *  0)]"'  Y  exp  [-uS(7)/c  j  (5a) 

v 

~r2 

T  3  I  °  [P(r  0]  1/Y  exp  (dS(7)/c  )  2r  '  dr  '  (5b) 

#0  v 

where  7*  is  the  Lagrangian  coordinate  of  a  particular  fluid  element, 
i.e.,  the  location  at  t  =  0  of  the  fluid  eler.er.t  that  is  at  r"  at  T  =  «, 
and  AS(?r)  is  the  increase  in  entropy  iensity  due  to  shocks  that  have 
traversed  this  fluid  element.  There  is  no  general  formula  for  1S(?'  , 
hut  in  Pef.  6  we  presented  heuristic  formulas  for  aSC?)  that  accurately 
described  the  outgoing  shock  ir.  the  case  of  a  Bennett  profile,  and  in 
fact  are  reasonably  accurate  for  the  first  outgoing  shock  in  the  general 
case.  F:r  the  ber.r.e*  *  profile,  If  3  0  for  the  centra-  channel,  since 
there  is  only  a  single,  outgoing  shock,  which  originates  at  7  *  0.8; 
thus  the  central  channel  is  formed  strictly  by  adiabatic  expansion,  and 
tr.e  main  effect  of  ’he  she  ex  is  to  broader.  *he  channel  wings.  Further¬ 
more,  the  time  evolution  is  very  smooth  ar.a  simple  for  the  Bennett  pro¬ 
file,  with  the  density  decreasing  nonotonically  tc  about  10X  belcw  the 
value  given  by  Eq.  5),  and  then  returning  slowly  tc  tha*  value. 

For  profiles  of  squarer  shape,  ve  have  seer,  that  the  hydrodynamic 
evolution  is  more  complicated,  i.e.:  i  A  strong  rarefaction  reduces 

the  density  to  a  value  well  below  that  of  Eq.(5),for  a  brief  period  vten 
~  The  temperature  is  also  at  a  minimum  at  this  tine,  (ii)  Beth 

the  iensity  ar.i  temperature  abruptly  increase  vr.er.  an  invar  i  shock  hits 
the  channel  axis.  This  shock  car.  be  quite  strong.  ill  There  is  a 
c*cd  leal  of  subsequent  structure  ir.  the  time  evolution.  Ir.  particular, 
the  incoming  shock  refle'ts  off  the  axis,  and  for  str-or.g  overpressures, 
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a  second  ingoing  shock  {which  is  always  weak)  is  observed,  (iv)  The 
channel  density  at  late  tines  when  pressure  balance  is  achieved  is  given 
by  Fq.  ;5)  and  is  significantly  decreased  below  the  adiabatic  value  by 
the  entropy  increase  uS  due  to  the  first  ingoing  shock,  e.g., 

■  D,  *t  *  •  )  *  O.Ohk  for  squarish  profiles  with  P  '  Tt  «s  compared 

o 

♦  an  adiabatic  value,  0.053.  The  channel  temperature  at  late  tine 
ieter-.ir.ed  by  pressure  balance'  is  always  much  larger  than  the  tempera- 


•  .re  luring  the  first  density  minimum,  e.g.  ^  =  63  at 


while 


.  ■  .  1  at  the  time  of  the  first  density  minimum. 

It  is  striking  that  these  effects  do  not  necessarily  require  that 
the  initial  overpressure  profiles  have  a  sharp  "comer"  or  discontinuity, 
but  ccur  also  for  profiles  that  are  only  moderately  squarer  ir.  shape 
than  the  Ser.nett  profile,  e.g.  the  generalised  Bennett  with  n  ■  -.5  (see 
Fig.  1  or  the  generalised  Gaussian  with  n  *  3  see  Fig.  ?).  In  fact 
the  ir.ward  shock  is  stronger  for  the  generalised  Gaussian  with  n  *  10 
than  for  squarer  profiles  with  n  >  ?C.  The  strong  rarefaction  and  in¬ 
var*.  shock  seem  to  be  due  primarily  tc  the  presence  of  a  relatively  flat 
-er.tral  regi  n  f  *he  profile.  This  can  be  understood  as  follows.  In 
a  ’uniform,  medium,  a  pressure  pulse  steepens  to  a  shock  as  a  result  of 
nonlinear  effects:  the  temperature,  and  therefore  the  local  sound  speed, 
is  larger  within  the  pulse,  and  therefore  the  center  cf  the  pulse  tends 
tc  overtak*  the  frrr.t.  This  tendency  is  reir.fcrcei  if  the  pressure  pulse 
*  -  propagat,  ing  down  a  temperature  gradient,  but  is  weakened  if  the  pulse 
-.st  prcpagate  t-vard  a  regicr.  of  increasing  temperature;  ir.  the  latter 
■ase,  the  front  of  the  disturbance  tends  to  r;r.  away  fror.  disturbance. 


As  a  result,  waves  rarefactions  or  compressions  propagating  up  a 
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temperature  gradient  tend  to  be  gentle,  and  compressions  tend  not  to 
steepen  Into  shocks.  Smooth  radial  pressure  profiles  such  as  the 
Bennett  profile  have  a  non-zero  Initial  temperature  gradient  except  on 
axis,  and  a  temperature  gradient  persists  during  the  Initial  expansion 
of  the  heated  region  (see  Fig.  8).  Thus  one  expects  subsequent  Inward- 
propagating  features  to  be  slow  and  gentle,  as  Is  Indeed  the  case:  the 
central  density  decreases  smoothly,  over  a  very  long  time  scale,  to  a 
value  slightly  below  Its  final  value,  and  returns,  even  more  slowly, to 
pressure  balance.  On  the  other  hand,  the  modified  Bennett  and  Gaussian 
distributions  with  n  >  5  have  a  pronounced  flat  central  region  (but  no 
real  corner  until  n  Is  much  larger),  which  persists  and  Is  conducive  to 
strong  Inward  propagating  rarefactions  and  shocks  at  early  times  (see 
Fig.  9).  For  very  square  profiles  (generalized  Gaussian  with  n  >10), the 
shock  strength  decreases  slightly.  On  the  other  hand,  an  initial  over¬ 
pressure  profile  In  the  form  of  a  truncated  Bennett  profile  does 
result  In  an  ingoing  shock  and  other  features  typical  of  a  "square” 
profile.  In  Fig.  10,  we  exhibit  the  evolution  resulting  from  an 
Initial  heating  profile 

1  ♦  P  (1  ♦  T2)-2  •  r  ^  ~ 

P  (?.r  -  0)  -  °  0  , 

1  •  r  ro 

wl£h  "  '4,  and  7  -  0.8.  Such  a  profile  does  not  have  an  Initially 
flat  central  region,  but  does  have  a  discontinuity  out  on  the  radial 
wings.  However,  we  see  In  Fig.  11  that  a  fist  central  temperature  pro¬ 
file  does  evolve  at  a  later  time,  due  to  the  rarefaction  wave  running 
In  from  the  discontinuity,  and  that  this  is  followed  by  the  formation 


of  an  ingoing  shock.  Thus  the  type  of  time  development  discussed  in 
this  paper  results  from  an  initial  overpressure  profile  that  is 
characterized  by  either  a  flat  central  region  or  a  discontinuity  at 
any  point. 

A.  Summary 

We  have  concluded  that  if  a  gas  is  heated,  by  a  directed  energy 
source,  instantaneously  on  a  hydrodynamic  time  scale,  in  such  a  way 
that  the  initial  overpressure  ptofile  ?(r,T)  "  0)  is  either  flat- 
topped  or  sharp-edged,  then  the  subsequent  expansion  will  not  be 
smooth  and  adiabatic  in  the  central  region,  as  it  is  for  a  Bennett 
pressure  profile. 6  Rather,  there  will  be  a  strong  rarefaction, 
followed  by  an  Inward  propagating  shock,  and  subsequent  oscillation 
or  "ringing"  of  the  channel  will  occur.  There  is  a  transient  period. 
Just  before  the  arrival  of  the  first  Inward  shock,  when  the  central 
ch*..uel  is  at  low  temperature  and  low  density.  For  square  deposition 
profiles,  the  channel  density  at  late  time,  when  pressure  balance 
with  the  surrounding  gas  is  achieved,  is  decreased  by  the  shock 
heating  to  about  80X  of  what  it  would  be  for  adiabatic  channel 
formation.  The  strength  of  the  shock  depends  mainly  on  the  shape 
of  the  profile.  The  time  scale  for  the  rarefaction  and  shock  depends 
mainly  on  the  strength  of  the  overpressure.  These  dependences  are 
extensively  illustrated  in  the  numerical  solutions. 
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Fig.  4  —  The  time  evolution  of  p(t~  •  0)  for  generalized  GauMian  heating  profiles 
with  overpressure  PQ  »  37.  and  for  various  values  of  the  squareness  index  n.  For 
clarity,  the  horizontal  axis  is  displaced  for  the  various  curves. 


Fig.  6  —  The  time  evolution  of  o{r~  •  0)  for  generalized  Bennett  heating  profiles 
with  overpressure  PQ  •  74.  and  for  various  values  of  the  squareness  index  n.  For 
clarity,  the  horizontal  axis  is  displaced  for  the  various  curves. 


